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ABSTRACT 

Classification theorems for linear differential equations in two real 
variables, possessing eigenf unctions in the form of the polynomi- 
als (the generalized Bochner problem) are given. The main result 
is based on the consideration of the eigenvalue problem for a poly- 
nomial elements of the universal enveloping algebras of the alge- 
bras sZ 3 (R), sZ 2 (R) © sZ 2 (R) and gl 2 (R) x R r+1 ,r > taken 
in the "projectivized" representations (in differential operators 
of the first order in two real variables) possessing an invariant 
subspace. General insight to the problem of a description of lin- 
ear differential operators possessing an invariant sub-space with 
a basis in polynomials is presented. Connection to the recently- 
discovered quasi-exactly-solvable problems is discussed. 
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Take the eigenvalue problem 

Ttp(x,y) = eip(x,y) (0) 

where T is a linear differential operator of two real variables, x, y and e is the 
spectral parameter. 

Definition. Let us give the name of the generalized Bochner problem 
to the problem of classification of the differential equations (0) possessing 
several eigenfunctions in the form of polynomials. 

In the paper a general method has been formulated for generating 
eigenvalue problems for linear differential operators, linear matrix differential 
operators and linear finite-difference operators in one and several variables 
possessing polynomial solutions. The method was based on considering the 
eigenvalue problem for the representation of a polynomial element of the 
universal enveloping algebra of the Lie algebra in a finite-dimensional, 'pro- 
jectivized' representation of this Lie algebra [|T|. 

In two previous papers 0, |3| it has been proven that in this approach 
the consideration of algebras s^R), s/2(R) q , osp(2, 2) in projectivized rep- 
resentations (in differential operators of the first order) provides in general 
both necessary and sufficient conditions for the existance of polynomial eigen- 
functions for ordinary linear differential operators, finite-difference operators 
in one variable and differential operators in two variables (one real and 
one Grassmann) (or, equivalently, 2x2 matrix differential operators) ]3[], 
respectevely. Particularly, it manifested the classification theorems, which 
imply the solution of the original Bochner problem (1929) posed for ordinary 
differential equations. In the present paper similar classfication theorems will 
be given for finite-order linear differential operators in two real variables in 
connection to the algebras s/ 3 (R), s/ 2 (R) © sZ 2 (R) and g/ 2 (R) X R m • Also 
an outlook on the problem for a general linear differential operators will be 
given. 

For future considerations, we define the space of all polynomials in x, y 
of finite degree as 

V N , M = (l,x 1 ,y 1 ,x 2 ,xy,y 2 ,...,x N y M ) (1) 

where N, M are non-negative integers, x, y G R are real variables. It is 
convenient to illustrate the space of polynomials of finite degree through the 
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Newton diagrams. In order to do this, let us consider two-dimensional plane 
and each point with the integer coordinates (k,p) put in correspondance with 
the monomial x k y p . 



1 Polynomials in two real variables 
1.1 Polynomials of the first type 

Now let us describe the projectivized representation of the algebra s^R) 
in the differential operators of the first order acting on functions of two real 
variables. It is easy to show that the generators have the form [|J 

Jl = y 2 d y + xyd x - ny , j\ = x 2 d x + xyd y - nx , 

■J3 = ydx j J\ = 9 X , Ji = dy , J2 — xdy , 

Jd = yd y + 2xd x - n , J d = 2yd y + xd x - n , (2) 

where x, y are the real variables and n is a real number. If n is a non-negative 
integer, the representation becomes finite-dimensional of the dimension (1 + 
n)(l + n/2). The invariant sub-space has a polynomial basis and is presented 
as a space of all polynomials of the following type 

P<P = (1; x, y; x\ xy, y 2 ;...; x n , x^y, . . . , xy n ~\ y n ) (3) 

or, graphically, (3) is given by the Newton diagram of the figure 1. 




Fig. 1 . Graphical representation (Newton diagram) of the space 

(see (3)). 
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Definition. Let us name a linear differential operator of the k-th order 
a quasi-exactly-solvable of the first type,T| (x, y), if it preserves the 
space V^p ■ Correspondingly, the operator E^. (x, y) G T^\x, y), which 
preserves the infinite flag V ( I] C V? C#C...C V® C ... of spaces of 
all polynomials of the type (3) , is named an exactly-solvable of the first 
type. 

LEMMA 1.1 Take the space T>jp . 

(i) Suppose n > (k — 1). Any quasi-exactly-solvable operator of the first 
type of the k-th order t£ (x, y), can be represented by a k-th degree polyno- 
mial of the generators (2). Ifn < (k—1), the part of the quasi-exactly-solvable 
operator T^\x,y) of the first type containing derivatives in x,y up to the 
order n can be represented by a n-th degree polynomial in the generators (2). 

(ii) Inversely, any polynomial in (2) is a quasi-exactly solvable operator 
of the first type. 

(Hi) Among quasi- exactly-solvable operators of the first type there exist 
exactly-solvable operators of the first type E^\x,y) C T^\x,y). 

Comment 1. If we define the universal enveloping algebra U g of a Lie 
algebra g as the algebra of all polynomials in generators, then the meaning of 
the Lemma is the following: T^x, y) at k < n+l is simply an element of the 
universal enveloping algebra U s i 3 (r\ of the algebra 3/3 (R) in representation 
(2). If k > n + 1, then Tk(x, y) is represented as a polynomial of the n-th 
degree in (2) plus B dxn ^ +lQym , where m = 0, 1, ... (n + 1) and B is any 
linear differential operator of the order not higher than (k — n — 1). 

Proof. The proof is based on the proof of irreducibility (3) and then on 
the application of Burnside theorem Q . 

3 I appreciate to J.Frohlich and M.Shubin for suggestion this way of the proof, which 
is easy and natural unlike the technical proof originally given at [^J. Similarly, one can 
prove the analogous lemmas for the case of linear differential operators in one variable, 
(see Lemma 1 in Q), finite-difference operators in one variable (Lemma 4 in and 
differential operators in one real and one Grassmann variable (Lemma 1 in [||) 
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Let us introduce the grading of the generators (2) in the following way. 
The generators are characterized by the two-dimensional grading vectors 

deg(Jl) = (+1,0), deg(J^) = (0,+l) , 

deg{Jl) = (-1,+1) , deg{Jl) = (+1,-1) , 

deg(J d ) = (0,0) , deg(J d ) = (0,0) , 

deg(Jf) = (0, -1) , de<?( J 2 2 ) = (-1, 0). (4) 

Apparently, the grading vector of a monomial in the generators (2) can be 
defined by the grading vectors of the generators by the rule 

degT = de^[(J 3 1 ) ni3 (J 2 1 )" 12 (J|)^(J|) nM (J d )^(J (i )^(Jf)" 31 (J?) n21 ] = 

("13 + "32 - "23 - "31 , "12 + "23 - "32 - "21) = (deg x T , deg y T) (5) 
Here the n's can be arbitrary non-negative integers. 

Definition. Let us name the grading of a monomial T in generators (2) 
the number 

deg(T) = deg x T + deg y T . 

We will say, that a monomial T possesses positive grading, if this number is 
positive. If this number is zero, then a monomial has zero grading. The no- 
tion of grading allows one to classify the operators Tk(x, y) in a Lie-algebraic 
sense. 

LEMMA 1.2 A quasi-exactly-solvable operator C £4z 3 (r) either has 
no terms of positive grading, iff it is an exactly- solvable operator. 

It is worth noting that among exactly-solvable operators there exists a 
certain important class of degenerate operators, which preserve an infinite 
flag of spaces of all homogeneous polynomials 

V^ = (x n ,x n ~ 1 y,...,xy n -\y n ) (6) 

(represented by a horizontal line in Fig.l). 

LEMMA 1.3 Linear differential operator T k (x , y) preserves the infinite 
flag Vq 1 " 1 C C C . . . C V^p C . . . of spaces of all polynomials of the 
type (6), iff it is an exactly-solvable operator having terms of zero grading 
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only. Any operators of such a type can be represented as a polynomial in 
the generators Jf , Jf , Jd, Jd (see (2)), which form the algebra S03 (BR. If 
such an operator contains only terms with zero grading vectors, this operator 
preserves any space of polynomials. 

THEOREM 1.1 Let n be non-negative integer. In general, the eigen- 
value problem for a linear symmetric differential operator in two real variables 
T k (x,y): 

T k (x,y)<p(x,y) = e<p(x,y) (7) 

has (n + l)(n/2 + 1) eigenfunctions in the form of a polynomial in variables 
x,y belonging to the space (3), iff T^.(x, y) is a symmetric quasi- exactly- 
solvable operator of the first type. The problem (7) has an infinite sequence 
of eigenfunctions in the form of polynomials of the form (3), iff the operator 
is a symmetric exactly-solvable operator. 

This theorem gives a general classification of differential equations 

k i=m d m Lo(r v) 

m=0 i=0 UX U y 

having at least one eigenfunction in the form of polynomial in x, y of the type 
(3). In general, the coefficient functions a>i£l-i(x,y) have quite cumbersome 
functional structure and we do not display them here (below we will give 
their explicit form for T2(x,y)). They are polynomials in x,y of the order 
(k + m) and always contain a general inhomogeneous polynomial of the order 
m as a part. The explicit expressions for those polynomials are obtained by 
substituting (2) into a general, the k-th order polynomial element of the uni- 
versal enveloping algebra t/ s z 3 (R) of the algebra 5/3 (R) . Thus, the coefficients 
in the polynomials di^-i(x,y) can be expressed through the coefficients of 
the k-th order polynomial element of the universal enveloping algebra U s i 3 (r)- 
The number of free parameters of the polynomial solutions is defined by the 
number of parameters characterizing a general, k-th order polynomial ele- 
ment of the universal enveloping algebra £/ s z 3 (r). In counting free parameters 
a certain ordering of generators should be fixed to avoid double counting 
due to commutation relations. Some relations between generators should be 
taken into account, specifically for the given representation (2), like 

Jl^d 2J 2 J(2 S^t/^ = T1J2 
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Jp d + 4J d - tJlJl = (n + S)Jl 

Jpd + Jlh - ZJUl = (n + 3)Ji 
3(J 2 1 J 1 2 + JlJl + J 2 3 J 3 2 ) + J d J d +J d J d - J d J d = 3J d + 3n + n 2 
2J d J d + 2J d J d - hJ d J d + 9JfJ 3 2 = (n + Q)J d + (n-3)J d + n 2 + 3n 
JdJd + JdJd + 5(J 2 1 J 1 2 + Jijj 5 ) + 2J 2 3 J 3 2 = (n + 3)(J d - J d ) 

J d J 2 — 2J d J 2 — 3J 2 jf = nJ 2 

J d Jl - 2J d Jl - 3Jpl = nJl (9) 

between quadratic expressions in generators (and the ideals generated by 
them). For the case of exactly-solvable problems, the coefficient functions 
a ijm-i( x iy) take the form 

p+q<m 

ai$-i(x,y) = a i,m-i, P , q x p y q (10) 

p,q=0 

with arbitrary coefficients. 

Now let us proceed to the case of the second-order differential equations. 
The second-order polynomial in the generators (2) can be represented as such 

T 2 = c aPa5 jp] + c Q/3 J« + c , (11) 
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where we imply summation over all repeating indices; a, /3, 7, 5 correspond to 
the indices of operators in (2) and for the Cartan generators we suppose both 
indices simulate d or d, all c's are set to be real numbers. Taking (9) into 
account, it is easy to show, that T 2 is characterized by 36 free parameters. 
Substituting (2) into (11), we obtain the explicit form of the second-order- 
quasi-exactly-solvable operator 

TfW) = 



[x 2 P™(x,y) + xP%(x,y) + P%(x, y)]-^ + 



d 2 

[xyP^(x,y) + Ps y i(x,y) + ^0(^,2/)]^^ + 



[y 2 P%(x,y) + yP%(x,y) + Pftfay)]^ + 



[xP* 2 (x,y) + P^(x,y) + P* (x,y)]-^ + 



[yPl 2 (x,y) + P^fay) + Pl ( x ,y)]^-+ 



[Pl 2 (x,y) + P? tl (x,y) + Pgjo&y)] (12) 

where P£ m (x,y) and P£ m (x,y) are homogeneous and inhomogeneous poly- 
nomials of the order k, respectively, the index m numerates them, the super- 
script c characterizes the order of derivative, which this coefficient function 
corresponds to. For the case of the second-order-exactly-solvable operator 
E 2 \x,y), the structure of coefficient function is similar to (12), except for 
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the fact that all tildeless polynomials disappear. If we denote the number of 
free parameter of the operator T with the symbol par(T), then it is easy to 
show that Q 

par(T$p(x,y)) = 36 . 

while for the case of an exactly-solvable operator (an infinite sequence of 
polynomial eigenf unctions in (8)) f] 

par (E^ix, y)) — 25 . 

An important particular case is when the quasi-exactly-solvable operator 
T2 (x, y) possesses two invariant sub-spaces of the type (3). This situation 
is described by the following lemma: 

LEMMA 1.4 Suppose in (11) there are no terms of grading 2: 

Cl2,12 = Ci 3 ,13 = Ci2,13 = (13) 

and if there exists some coefficients c 's and a non-negative integer N such 
that the conditions 

C12 = (n - N - m)ci2,d + (n - 2N + m)c 12 j + (N - m)ci 2 ,32 , 

C13 = (n - N - m)c 13jd + (n - 2N + m)c 13 j + mci 3i2 3 , (14) 

are fulfilled at all m = 0, 1, 2, ... , N, then the operator T% (x, y) preserves 
both V n andVjq, and par (T2 T \x,y)) = 31. 

Now let us proceed to the important item: under what conditions on the 
coefficients in (11) the second-order-quasi-exactly-solvable operators can be 

4 Recall, that for the case of the second-order-quasi-exactly-solvable differential operator 
in one real variable, the number of free parameters was equal to 9 (see||), for the case of 
one real and one Grassmann variables this number was 25 (see ||). 

5 Recall, that for the case of the second-order-exactly-solvable differential operator in 
one real variable, the number of free parameter was equal to 6 (see Q), for the case of one 
real and one Grassman variables this number is 17 (see ||). 
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reduced to a form of the Schroedinger operator after some gauge transfor- 
mation 

f(x,y)e t ^T 2 (x,y)e- t ^ = —A g + V(x,y) (15) 

where /, t, V are some functions in R 2 , and A g is the Laplace-Beltrami oper- 
ator with some metric tensor g^; fi, v — 1,2 0. This is a difficult problem for 
which there is yet no complete solution. In Refs. j|, a few multi-parametric 
examples were constructed. However, there is a quite general situation, for 
which a rather wide class of the solutions can be obtained. 

The algebra 5/3 (R) in realization (2) contains the algebra sc»3(R) as a 
sub- algebra 

J 1 = (1 + y 2 )d y + xyd x — ny , J 2 = (1 + x 2 )d x + xyd y — nx , 

J 3 = xd y - yd x (16) 

If the parameter n is a non- negative integer (and coincides with that in (2)), 
the same finite-dimensional invariant sub-space (see (3) and Fig.l) as in 
the original s^R) occurs. For this corresponding finite-dimensional 

representation (3) is reducible and unitary. It has been proven @, 0], that any 
symmetric bilinear conbination of generators (16), T 2 = c a pJ a J^, c Q/ g = cp ai 
can be reduced to a form of the Laplace-Beltrami operator plus a scalar 
function. In general, the metric tensor g^ u is not degenerate f\. Generically, 
the potential V(x, y) is given by a rational function and has no dependence on 
the spin n (see (16)). However, if the matrix c Q/ g has one vanishing eigenvalue, 
a certain mysterious relation appears 0, U 

V(x,y,{c}) = ^-R(x,y,{c}), (17) 
16 

where R(x,y, {c}) is the scalar curvature calculated through the metric ten- 
sor attached in the Laplace-Beltrami operator. This can imply that the 
corresponding Schroedinger operator has a purely geometrical nature! The 
real meaning of this fact is still not understood. 

So, we arrive to the two-dimensional exactly-solvable Schroedinger equa- 
tions. As a consequence of the quadratic Casimir operator for sos(R) in 

6 Our further consideration will be restricted the case / = 1 only 
7 Degeneracy of occurs, for example, if c a p has two vanishing eigenvalues 
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the form (16) being non-trivial and commuting with T 2 , the functional space 
of T2 is subdivided into the finite-dimensional blocks corresponding to the 
irreducible reps of so 3 (R). 

1.2 Polynomials of the second type 

In Ref. [0] we studied the quasi-exactly-solvable operators in one real vari- 
able. It turned out that the solution to this problem was found using a 
connection with the projectivized representation of the algebra s^R)- As a 
natural step in developing the original idea, let us consider the projectivized 
representation of the direct sum of two algebras s/ 2 (R)- 

The algebra s/ 2 (R) ©s/ 2 (R) taken in projectivized representation acts on 
functions of two real variables. The generators have the form (see e.g.Q) 

= x 2 d x - nx , Jy = y 2 d y - my 

where x, y are the real variables and n, m are non-negative integers. There 
exists the finite-dimensional representation of dimension (n + l)(m + 1). Ev- 
idently, the invariant sub-space has a polynomial basis and is presented as 
a space of all polynomials with the Newton diagram shown in Fig. 2. We 

denote this space as Vi 1 ^. 

^ n > m x n y m 




Fig. 2. Graphical representation (Newton diagram) of the space V^m ■ 
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Definition. Let us name a linear differential operator of the k,p-ih order, 
containing derivatives in x and y up to k— th and p— th orders, respectively 0, 
a quasi-exactly-solvable of the second type, T^^ (x,y), if it preserves 

the space V^. Correspondingly, the operator E^y (x, G T^(x, y), 
which preserves either the infinite flag Pq 7 ^ ^ "Pi,™ c ^2^? C . . . C C 
. . . , or the infinite flag 7^ c V n,x C ^n,P C. . . C Pg, C ... of spaces 
of all polynomials, is named an exactly-solvable of the type 2 X or 2 y , 
respectively. . 

LEMMA 2.1 Take the space V}$. 

(i) Suppose n > (k — 1) and m > (p — 1). yln?/ quasi-exactly-solvable 
operator of the second type (k,p)-th order Tj^ (x, y), can be represented by 
a (k,p)-th degree polynomial of the generators (18). If n < (k — 1) and/or 
m < (p — 1), the part of the quasi- exactly-solvable operator T^p(x, y) of the 
second type containing derivatives in x, y up to the order n, m, respectively, 
can be represented by a (n,m)-th degree polynomial in the generators (18). 

(ii) Inversely, any polynomial in (18) is a quasi- exactly solvable operator 
of the second type. 

(Hi) Among quasi- exactly-solvable operators of the second type there exist 
exactly-solvable operators of the second type ~E^(x,y) C T^' (x,y). 

The proof is analogous to the proof of Lemma 1.1 and is based on the 
irreducibility of V^m an d Burnside theorem. 

Similarly, as has been done for the algebra s^R), one can introduce the 
notion of grading: 

deg(J+) = (+1,0) , deg(J+) = (0,+l) , 

deg(J° x ) = (0,0) , deg(J° y ) = (0,0) , 

deg(J-) = (-1, 0) , deg(J y ) = (0, -1). (19) 

8 So a leading derivative has a form ^p^p- ■ Also we will use a notation through this 
section T N (x, y) implying that in general all derivatives of the order TV are presented 
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Apparently, the grading vector of a monomial in the generators (18) can be 
defined by the grading vectors (19) of the generators by the rule 

degT = ^[(J+)-+(J°)"-(J-)--(J, + )^(jJ)^'(J;)"-] = 

{n x+ - n x _ , n y+ - n y _) = (deg x (T) , deg y {T)) (20) 

Here the n's can be arbitrary non-negative integers. 

Definition. Let us name the x— grading, y— grading and grading of a 
monomial T in generators (18) the numbers deg x (T), deg y (T) and deg{T) = 
deg x (T) + deg y (T), respectively We say that a monomial T possesses positive 
x— grading (y— grading, grading), if the number deg x (T) (deg y (T), deg(T)) 
is positive. If deg x (T)(deg y (T),deg(T)) = 0, then a monomial has zero 
x— grading (y— grading, grading). 

The notion of grading allows one to classify the operators (x, y) in a 
Lie-algebraic sense. 

LEMMA 2.2 The quasi-exactly-solvable operator T^t ,(x,y) of the sec- 
ond type preserves the infinite flag c ^i,m c ^2,2 C . . . C C . . . 
of spaces of all the polynomials, iff it is an exactly-solvable operator of the 
type 2 X having no terms of positive x-grading, deg x > 0. 

The quasi-exactly-solvable operator T^p(x,y) preserves the infinite flag 

V^q C V^P C V^P C ... C C ... of all spaces of the polynomials, iff 

it is an exactly-solvable operator of the type 2 y having no terms of positive 
y -grading, deg y > 0. 

// a quasi-exactly-solvable operator of the second type contains no terms 
of positive grading, this operator preserves the infinite flag Vq 1 ^ C C 
C . . . C Vj:p C . . . of spaces of all the polynomials of the type (3) and is 
attached to the exactly-solvable operator of the first type. 

THEOREM 2.1 Let n,m be non-negative integers. In general, the 
eigenvalue problem (7) for a linear symmetric differential operator in two 
real variables T k>p (x,y) has (n + l)(m + 1) eig en functions in the form of a 
polynomial in variables x, y belonging to the space V^J^ , iff T k:P (x, y) is a 
quasi-exactly-solvable symmetric operator of the second type. The problem 
(7) has an infinite sequence of eig en functions in the form of polynomials 
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belonging the space V^m at fixed m (n), iff the operator is an exactly- solvable 
symmetric operator of the type 2 X (2 y ). 

This theorem gives a general classification of differential equations (8), 
having at least one eigenfunction in the form of a polynomial in x, y of the 
type V^V. In general, the coefficient functions a^_j(x, y) in (8) have a quite 
cumbersome functional structure and we do not display them here (below we 
will give their explicit form for T 2 (x, y)). They are polynomials in x, y of 
the order (k+m). The explicit expressions for those polynomials are obtained 
by substituting (18) into a general, the k-th order polynomial element of 
the universal enveloping algebra £/sz 2 (r)© s j 2 (r) °f the algebra 3/2 (R) © ^(R) 
. Thus, the coefficients in the polynomials af^L^x, y) can be expressed 
through the coefficients of the fc-th order polynomial element of the universal 
enveloping algebra ?7 s i 2 (R)© s / 2 (R) • The number of free parameters of the 
polynomial solutions is defined by the number of parameters characterizing 
a general, fc-th order, polynomial element of the universal enveloping algebra 
U s i 2 (n)<$ s i 2 (R). In counting free parameters a certain ordering of generators 
should be fixed to avoid double counting due to commutation relations. Also 
some relations between generators should be taken into account, specifically 
for the given representation (18), like @ 

^x^x ~ Jx^x + {~2 ^)^x = ~ ~2 (~2 ~ ^ 

Jpy ~ JyJy + (y + W = " V ( 21 ) 

between quadratic expressions in generators (and the ideals generated by 
them) f\. 

Now let us proceed to the case of the second-order differential equations. 
The second-order polynomial in the generators (18) can be represented as 
such 

ji r xx T a T@ -I- r xy T a Jf 3 -I- r yy T a 1$ -I- r x T a -I- r y T a -I- r (99~\ 

±2 — L -al3 J x J x ' c a/3 J x J y ' °af3 J y J y ' c a J x ' °a J y ' c ) K**) 

where we imply summation over all repeating indices and a, (3 = ±, 0; and all 
c's are set to be real numbers. Taking (21) in account, it is easy to show that 

9 For this case they correspond to quadratic Casimir operators 
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T 2 is characterized by 26 free parameters. Substituting (18) into (22), we 
obtain the explicit form of the second-order-quasi-exactly-solvable operator 



T { 2 H \x,y) 



few + y p U x )]^ + \- p Uv) + xP Uv)\^y + p U^y) ( 23 ) 

where P^ m (x,y) and Pj: m (x,y) are homogeneous and inhomogeneous poly- 
nomials of the order k, respectively, the index m numerates them, and the 
superscript c characterizes the order of the derivative corresponding to this 
coefficient function. For the case of the second-order-exactly-solvable opera- 
tor of 2 X type 

Qv>W-j& + Wife) + QToiv)]^ + QfM^, + 

Ql^x) + yQto(x)}^ + QIM + xQ%{y)\^ + Q%(y) (24) 

For the case of the second-order-exactly-solvable operator of 2 y type the 
functional form is similar to (24) with the interchange x <-> y. 

It is easy to show that the number of free parameters are 

par(T { 2 n) (x,y)) = 26 . 



10 Recall, that for the case of the second-order-quasi-exactly-solvable differential operator 
in one real variable, the number of free parameters was equal to 9 (see[||), for the case 
of one real and one Grassmann variables this number was 25 (see ||). Also, for the case 
of the second-order-exactly-solvable differential operator in one real variable, the number 
of free parameter was equal to 6 (see 0), for the case of one real and one Grassmann 
variables this number is 17 (see |||) (see also p. 9 for comparison to the case of the first 
type polynomials). 
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par(E ( 2 II] (x,y)) = 20 . 

Similar to the previous cases, there is a very important particular case of 
quasi-exactly-solvable operators of the second order, where they possess two 
invariant sub-spaces. 

LEMMA 2.3 Suppose in (22) there are no terms of x— grading 2: 

c x + x + = (25) 

and if there exists some coefficients c's and a non-negative integer N such 
that the conditions 

c x A = o , 

4- = , 

c\ = (n/2 - N)c% + (m/2 - fc)c& , (26) 

are fulfilled at all k = 0, 1, 2, ... , m, then the operator T 2 I \x, y) preserves 
both V^m and V}J m , besides that par(T 2 ( x i V)) = 22. 

Generically, the question of the reduction of quasi-exactly-solvable oper- 
ator of the second type T 2 T \x, y) to the form of the Schroedinger operator 
is still open. In the papers [|], || several multi-parametrical families of those 
Schroedinger operators were constructed. As in the case of the first type 
quasi-exactly-solvable operators, corresponding Schroedinger operators con- 
tain in general the non-trivial Laplace-Beltrami operator. 

The above analysis of linear differential operators preserving the space 
V^m can be naturally extended to the case of linear finite-difference operators 
defined through the Jackson symbol D: 

n ,, x f{x) - f{qx) . , 

Df(x) = -— r + f(qx)D 

(I — q)x 

where f(x) is a real function and q is a number, instead of continuous deriva- 
tive. All above-described results hold with replacement of the algebra 
s/ 2 (R) © sl 2 (R) to the quantum algebra s/ 2 (R-) g © sl 2 (R) q in 'projectivized' 
representation (see|jl], @). In the limit q — > 1 all results, which can be ob- 
tained, coincide to the results of present Section. 

11 with minor modifications 
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1.3 Polynomials of the third type 



The third case, which we are going to discuss here, corresponds to the Lie 
algebra g^CR) X R r+1 (semidirect sum of ^(R) with a (r + 1)— dimensional 
abelian ideal; the Case 24 in the classification given the paper ||). This 
family of the Lie algebras, depending on an integer r > 0, can be realized in 
terms of the first order differential operators 



J 1 = d. 



X ) 



J 2 = xd x - ^ , J 3 = yd y - — , 
3 3r 

J 4 = x 2 d x + rxydy — nx , 

j 5+i = x { d y , z = 0,1,..., r , (27) 

where x, y are the real variables and n is a real number q If n is a non- 
negative integer, the representation becomes finite-dimensional. The invari- 
ant sub-space has a polynomial basis and is presented as a space of all poly- 
nomials of the form 

i+rj<n 

-P% 1} = E <H**V (28) 

i,j>0 

or, graphically, (28) is given by the Newton diagram of the figure 3. The 
general formulas for the dimension of corresponding finite-dimensional rep- 
resentation (28) is given by 

dimV irn) = ^ + (T +£ n + ^1 (29) 



where for small r 



012;, 



«3,n 



OCl, n = * , 

J 3 at odd n 
I 4 at even n 

4 at (n+1) multiple 3 
6 at other n 



12 It is worth noting that at r = 1 the algebra {c^ 2 (R-) X R 2 } C sZ 3 (R). Thus, this 
case is reduced to one about the first type polynomials (see Ch.1.1). Hereafter, we include 
r = 1 into consideration just for the purpose of completeness. 
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5 at (n+1) multiple 4 
«4,n = { 8 at other n 

9 at (n+3) multiple 4 

Actually, the numbers a r ^ n are given by all possible products m(2r — m) at 
m — 1, 2, . . . , r. 




Fig. 3 . Graphical representation (Newton diagram) of the space V^ 1 ^ 

(see (28)). 

Definition. Let us name a linear differential operator of the k-th order a 
quasi-exactly-solvable of the r— third type, T ( [' III \x,y), if it preserves 
the space V^ 1 ^ ■ Correspondingly, the operator ~E ( £' III \x, y) G T^' IH \x,y), 

which preserves the infinite flag V^o^ C V^l^ C V^ 1 ^ C . . . C V^ 1 ^ C . . . 
of spaces of all polynomials of the type (28) , is named an exactly-solvable 
of the r— third type. 

LEMMA 3.1 Take the space V^ J) . 

(i) Suppose n > (k—1). Any quasi- exactly-solvable operator of the r— third 
type of the k-th order T^' IH \x,y), can be represented by a k-th degree poly- 
nomial of the generators (27). If n < (k — 1), the part of the quasi- exactly- 
solvable operator T ( [' IH \x,y) of the r— third type containing derivatives in 
x, y up to the order n can be represented by a n-th degree polynomial in the 
generators (27). 

(ii) Inversely, any polynomial in (27) is a quasi-exactly solvable operator 
of the r— third type. 
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(Hi) Among quasi-exactly-solvable operators of the r— third type there exist 
exactly- solvable operators of the r— third type ~Ei k r,III \x,y) C T k r ' IIT \x, y). 

The proof is analogous to the proof of the Lemma 1.1 and 2.1. It is based 
on irreducibility V^ 1 ^ and Burnside theorem. 

One can introduce the grading of the generators (27) in analogous way 
as has been done before for the cases of the algebra s/s(R) (see (4)) and 
s/ 2 (R) © s^(R) (see (19)). All generators are characterized by the two- 
dimensional grading vectors 

de^(J 1 ) = (-l,0) , 

deg(J 2 ) = (0,0) , deg(J 3 ) = (0,0) , 
deg(J 4 ) = (1,0) , 

deg(J 5 ) = (0, -1) , deg(J 6 ) = (1, -1) , . . . , deg(J 5+r ) = (r, -1) (30) 

Similarly as before, the grading vector of a monomial in the generators (27) 
can be defined through the grading vectors of the generators (27) (cf.(5),(20)). 

Definition. Let us name the grading of a monomial T in generators (27) 
the number 

deg(T) = deg x T + rdeg y T 

(cf. the case of sZ 3 (R)). We will say that a monomial T possesses positive 
grading if this number is positive. If this number is zero, then a monomial 
has zero grading. The notion of grading allows one to classify the operators 
T k (x, y) in a Lie-algebraic sense. 

LEMMA 3.2 A quasi- exactly- solvable operator T^ 111 ^ C U g i 2 ^ jxr^+i 
either has no terms of positive grading, iff it is an exactly-solvable operator. 

THEOREM 3.1 Let n and (r — 1) be non-negative integers. In gen- 
eral, the eigenvalue problem (7) for a linear symmetric differential operator 
of the k—th order in two real variables Tk{x,y) has a certain amount of 
eigenf unctions in the form of a polynomial in variables x, y belonging to the 
space (28), iffT k (x,y) is a quasi-exactly-solvable, symmetric operator of the 
r— third type. The problem (7) has an infinite sequence of eigenf unctions 
in the form of polynomials belonging to (28), iff the operator is an exactly- 
solvable, symmetric operator of the r— third type. 
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This theorem gives a general classification of differential equations (8), 
having at least one eigenfunction in the form of a polynomial in x, y of the 
type V^n^ ■ m general, the coefficient functions af^^x, y) in (8) are polyno- 
mials in x, y and have a quite cumbersome functional structure and we do not 
display them here (below we will give their explicit form for T2' IH \x,y)). 
The explicit expressions for those polynomials are obtained by substitut- 
ing (27) into a general, k-th order polynomial element of the universal en- 
veloping algebra U g i 2 m\ c<Fr+i of the algebra g^CR) >< R r+1 . Thus, the 
coefficients in the polynomials af^_ { (x, y) can be expressed through the co- 
efficients of the /c-th order polynomial element of the universal enveloping 
algebra £^ 2 (r) xrh-i. The number of free parameters of the polynomial so- 
lutions is defined by the number of parameters characterizing general /c-th 
order polynomial element of the universal enveloping algebra. In counting 
free parameters a certain ordering of generators should be fixed to avoid 
double counting due to commutation relations. Also some relations between 
generators should be taken into account, specifically for a given representa- 
tion (27), like 

jV 5 - J 1 J 6 + ^J 5 = , 
J 1 J 4 - J 2 J 2 - rJ 2 J 3 - J 2 - r(| + l)J 3 = + , (31) 

71 

J 2 J Q+l + rJ 3 J 6+l - J 4 J 5+l ~ G + l)J 6+i = , 

o 

ati = 0, 1,2, ...,(r- 1) , (32) 

jljW _ J2j 6+i - (- + 1) J 6+l = , 

3 

ati = 0, 1,2, ...,(r-2) , (33) 

j5j7+i _ _ jb+k j7+i-k 

at k = 0, 1, 2, . . . , and 2k < (2 + i) and % = 0, 1, 2, . . . , (2r - 2) (34) 

between quadratic expressions in generators (and the ideals generated by 
them). 
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Now let us proceed to the case of the second-order differential equations. 
The second-order polynomial in the generators (27) can be represented as 

T 2 = c a pJ a jP + c a J a + c , (35) 

where a, (3 = 1, 2, . . . , (5 + r) and we imply summation over all repeating 
indices; all c's are set to be real numbers. Taking the relations (31)-(34) in 
account, one can count the number of free parameters and obtain 

par(Tp IH \x,y)) = 5(r + 4) (36) 

Substituting (27) into (35), we obtain the explicit form of the second-order- 
quasi-exactly-solvable operator 

Ti' HI \x,y) = 

[P&W ]§; + + V??M\§j + P U X ) ( 37 ) 

where P% (x, y) are inhomogeneous polynomials of the order k, the index m 
numerates them, and the superscript c characterizes the order of the deriva- 
tive corresponding to this coefficient function. For the case of the second- 
order-exactly-solvable operator of r-third type 

Ep ni \x,y) = 

P % {X) & + [P ^ {X) + yP ^ x)] ^dy- + [P1 ^ X) + yP ^ X ^ + 

In this case the number of free parameters is equal to 

par(Ep in \x,y)) = 5(r + 3) (39) 

(cf.(36)). 
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Similar to the previous cases, there is an important particular case of 
r-third type quasi-exactly-solvable operators of the second order, where they 
possess two invariant sub-spaces. 

LEMMA 3.3 Suppose in (35) there are no terms of grading 2: 

c 44 = (40) 

and if there exists some coefficients c's and a non-negative integer N such 
that the conditions 

N N 
Ca = m)c 2 4 + (y - k)c M , (41) 

are fulfilled at all m, k = 0, 1, 2, . . . such that m + rk = N , then the operator 
I2 {x, y) preserves bothV^J andVl 1 ^, and par(T2' IH \x,y)) = 5r + 17. 

Generically, the question of the reduction of quasi-exactly-solvable op- 
erator of the r-third type T^' 11 (x, y) to the form of the Schroedinger op- 
erator is still open. Initially || several multi-parametrical families of those 
Schroedinger operators were constructed. As in the case quasi-exactly-solvable 
operators of the first and the second type, corresponding Schroedinger oper- 
ators contain in general a non-trivial Laplace-Beltrami operator. 



1.4 Discussion 

Through all the above analysis, very crucial requirement played a role: we 
considered the problem in general position, assuming the spaces of all poly- 
nomials contain polynomials with all possible coefficients. If this requirement 
is not fullfilled, than the above connection of the finite-dimensional space of 
all polynomials and the finite-dimensional representation of the Lie algebra 
is lost and degenerate cases occur. This demands a separate investigation. 

So, we described three types of finite-dimensional spaces (see Fig. 1-3) 
with a basis in polynomials of two real variables, which can be preserved by 
linear differential operators. Very natural question emerges: is it possible to 
find linear differential operators, which preserve the space of polynomials of 
finite order other than shown in Fig. 1,2, 3 ? 
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Conjecture 1. If a linear differential operator Tk(x,y) does preserve 
a finite- dimensional space of all polynomials other than V„, V^ m , V^J , then 
this operator preserves also a certain infinite flag of finite- dimensional spaces 
of all polynomials. 

2 General consideration 

The main conclusion of the considerations of the present paper and the papers 
|I|, |2], |3| one can formulate as the following theorem: 

(Main) Theorem. Take a Lie algebra g, realized in differential opera- 
tors of the first order, which possesses a finite- dimensional irreducible repre- 
sentation V ( a finite module of smooth functions). Any linear differential 
operator, having the invariant subspace V, which coincides to the above finite 
module of smooth functions, can be represented by a polynomial in generators 
of the algebra g plus an operator annuhilating V . 

Proof. Use of the Burnside theorem. 

As a probable extension of this theorem, we assume that the following 
conjecture hold 

Conjecture 2. If a linear differential operator T acting on functions in 
R k does possess the only finite- dimensional invariant subspace with polyno- 
mial basis, this finite- dimensional space coincides to a certain finite- dimensional 
representation of some Lie algebra, realized by differential operators of the 
first order. 

The interesting question is how to describe finite modules of smooth func- 
tionsin R k , which can serve as invariant sub-spaces of linear operators. Ev- 
idently, if we could found some realization of a Lie algebra g in differential 
operators, possessing an irreducible finite module of smooth functions, one 
can immediately consider the direct sum of several species of g acting on 
the functions given at the corresponding direct products of spaces, where 
the original algebra g acts in each of them. As an example, this procedure 
is presented in the Section 1.2 : the direct sum of two s£ 2 (R) leads to the 
operators acting on the functions at R 2 and a rectangular VjJQ is preserved. 
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Taking the direct sum of k species of SI2 (R) acting at R k , we arrive to 
k- dimensional parallelopiped as the invariant sub-space. In this case the 
algebra (R) plays the role of a primary algebra giving rise to a multidi- 
mensional geometrical figure (Newton diagram) as the invariant sub-space 
of some linear differential operators. In the space R 1 there is only one such 
"irreducible" Newton diagram - a finite interval, in other words, a space of 
all polynomials in x of degree not higher than a certain integer. Thus, in R 2 
the rectangular Newton diagram is "reducible" stemming from the product 
of two intervals as one-dimensional Newton diagrams.. In R 2 the irreducible 
Newton diagrams are exhausted by different triangles, connected to the al- 
gebras 3/3 (R) and {g/2 (R) X R r+1 }. We could not find any other Newton 
diagrams in R 2 (see conjecture 1). 

Before to present our present knowledge about Newton diagrams for 
the R k case, let us recall on the regular representation (in the first order 
differential operators) of the Lie algebra sIn(R) given on the flag man- 
ifold, which acts on the smooth functions of N(N — l)/2 real variables 
z i,i+q ,i — 1) 2, . . • , (N — 1), q = 1, 2, . . . , (N — i). The explicit formulas 
for generators are given by (see fToH ) 

d N d 
D ^ = ^ E z i+hv 



u ^i,i+l q=i+2 u4, i,q 



i q i-1 g 

D(fi) = Z Q,i+l Z iA+l Tj I" X/( Z <7,i+l ~~ z q,i z i,i+l 



q=l ^^q,t+l q= \ ^^q,t 

/ j z i,q o ^i z i,i+l 
q=i+2 C>Z i+l,q 

N Q i-1 Q i 

D(hi) = - + E 757- ~ E n *-p ( 42 ) 

q=i+l UZ t,q q=l uz q,t p= \ 

where we use a notation e^, f\, hi for generators of positive and negative roots, 
and the Cartan generators, correspondingly. If rij are non-negative integers, 
the finite-dimensional irreducible representation of sZjv(R) will occur in the 
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form of inhomogeneous polynomials in variables Zi^ +q . The highest weight 
vector is characterized by the integer numbers ni, i — 1,2, ... , (N — 1). 

Connected to the space R 3 , at least, two irreducible Newton diagrams 
appear : a tetrahedron, connected to a degenerate finite-dimensional rep- 
resentation of sl±(R,) (highest weight vector has two vanishing integers and 
one non- vanishing) , and a certain geometrical figure, corresponding to the 
regular representation of 5/3 (R) given on the flag manifold. We do not know 
if those exhaust all irreducible Newton diagrams or none. For the general 
case, R k there also exist polyhedra connected to sZfc +1 (R) (except one, all 
integers, characterizing the highest weight vector are vanishing) and geo- 
metrical figures related to some finite-dimensional regular representations of 
sZfc + i_i(R), i > 0. It is easy to show since the algebra so fc+1 (R) C s/ fc+1 (R) 
and has the polyhedron as the invariant sub-space 0, any quadratic polyno- 
mial in generators of the algebra so^ + i(R) with real symmetric coefficients 
can be reduced to the form the Laplace-Beltrami operator plus a scalar func- 
tion by means of a gauge transformation (see discussion in Ref. [|7|]). 

The above procedure can be extended to the case of the Lie super-algebras 
and also the quantum algebras. For the former, for linear operators acting on 
functions of real and Grassmann variables one can appear finite-dimensional 
invariant subspaces others than finite-dimensional representations (see [§]). 
However, those subspaces are connected to the representations of polyno- 
mial elements of the universal enveloping algebra and, correspondingly, they 
can be constructed through finite-dimensional representations of the original 
super- algebra. For the latter, we could not find any irreducible (non-trivial) 
Newton diagrams in R k for k > 1; they will occur, if a quantum space is 
considered instead of ordinary one 0. Also one can construct quasi-exactly- 
solvable and exactly-solvable operators considering "mixed" algebras: a di- 
rect sum of the Lie algebras, Lie superalgebras and the quantum algebras, 
realized in the first order differential and/or difference operators, possessing 
a finite-dimensional invariant subspace. 

We have described above the linear operators, which possess a finite- 
dimensional invariant sub-space with a polynomial basis. Certainly, by means 

13 Corresponding finite-dimensional representation is unitary and reducible (see e.g. the 
discussion in the end of Section 1.1) 

14 I am grateful to O.V.Ogievetski for discussion of this topic. 
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of diffeomorphism and gauge transformation one can obtain the operators 
with an invariant subspace with a non-polynomial basis (but emerging from 
polynomial one). The open question is : is it possible to find the operators 
possessing a finite-dimensional invariant sub-space with some explicit basis, 
which can not be reduced to a polynomial one using a diffeomorphism and 
gauge transformation? Up to now no examples of this type have been found 
(see Conjecture 2). Since the original aim of this whole investigation was 



the construction of the quasi-exactly-solvable Schroedinger operators [11 



S H i) Hi; a resolution of this problem plays a crucial role in leading to a 
complete classification of the quasi-exactly-solvable Schroedinger operators. 
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